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Thermoelectric coefficients of a two dimensional electron gas (2DEG) with the Rashba spin-orbit 
interaction (SOI) are presented here. In absence of magnetic field, thermoelectric coefficients are 
enhanced due to the Rashba SOI. In presence of magnetic field, the thermoelectric coefficients of 
spin-up and spin-down electrons oscillate with different frequency and produces beating patterns 
in the components of the total thermoelectric power and the total thermal conductivity. We also 
provide analytical expressions of the thermoelectric coefficients to explain the beating pattern for- 
mation. We obtain a simple relation which determines the Rashba SOI strength if the magnetic 
fields corresponding to any two successive beat nodes are known from the experiment. 
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I. INTRODUCTION 

There has been a rapid growth of research interest on 
the Rashba SOI in low-dimensional condensed matter 
system after the proposal of spin field effect transistor by 
Datta and Das^. This is due to the possible applications 
in spintronics devices^r— . The SOI is responsible for other 
interesting effects like spin Hall effect^, spin dynamics 
and zitterbewegung^"? . In the narrow gap semiconductor 
heterostructures, the dominant Rashba SOI^'^*^ appears 
due to the asymmetric quantum wells. The Rashba SOI 
strength is proportional to the internally generated crys- 
tal field. This strength can also be enhanced by applying 
suitable electric field perpendicular to the plane of the 
electron's motio n^^i^^ . 

A pseudo Zeeman effect occurs at finite momentum of 
the electron due to the Rashba SOI even in absence of 
magnetic field. A direct manifestation of pseudo Zee- 
man effect due to the SOI is a regular beating pattern 
in the magnetoelectric transport measurements such as 
Shubnikov-de Hass (SdH) oscillationsi^ in 2DEG. These 
oscillations occur due to two closely spaced different fre- 
quency of spin-up and spin-down electrons. The Rashba 
SOI strength is determined by analyzing the beating pat- 
terns in the SdH oscillation a^'^i^^ . The SOI was deter- 
mined by fitting the experimental data with the model 
calculations for the SdH oscillations. Later, many realis- 
tic approach was considered and the estimated strength 
is in good agreement with the extrapolated resultsi^"— . 
Recently, there is an interesting proposal^ of determin- 
ing the Rashba SOI strength by analyzing the beating 
patterns in the Weiss oscillations^'^'^"'^. 

On the other hand, thermoelectric properties of 
materials^ have attracted considerable interest from 
both experimental and theoretical point of view due 
to potential applications in technolog y^ - t^^ . There is 
a strong effect of perpendicular magnetic field on ther- 
mal transport properties of any system. Therefore, the 
magnetothermal coefficients can be used as an additional 
probe. In presence of perpendicular magnetic field, the 
diffusing charge carriers experience the Lorentz force. 



This produces a transverse electric field in addition to the 
longitudinal electric field. The longitudinal thermopower 
or the Seebeck coefficient is defined as Sxx = 
On the other hand, the transverse thermopower or the 
Nernst coefficient is defined as Sxy = — Here, V14 
and VVjy are the induced voltage generated by the ther- 
mal gradient and the magnetic field, respectively. The- 
oretical and experimental studies on thermoelectric co- 
efficients of 2DEG systems in presence of magnetic field 
started after the discovery of the quantum Hall effect. In 
most of the thermoelectric measurements of 2DEG sys- 
tems, the thermopower is being measured since the ther- 
mal resistivity of a 2DEG is extremely high. The Nernst 
coefficient is quite sensitive to various properties of the 
systems e.g. shape of the Fermi surface as well as electron 
mean free patb^. It is being used as a probe to study var- 
ious strongly correlated electron systems such as Kondo 
lattices^^ and graphene field effect transistors^^'^*. More- 
over, the thermopower S and the thermal conductivity k 
are used as the metrics to measure the thermoelectric 
performance^^. In addition to these, we will show here 
that magntethermoelectric coefficients can also be used 
to determine the Rashba SOI strength. 

There are mainly two mechanisms contribute to the 
thermal conductivity and the thermopower, namely the 
thermodiffusion and phonon drag. Generally, the phonon 
drag contribution is vanishingly small at very low temper- 
ature. In absence of the magnetic field, the diffusive ther- 
mopower has been continuously reported in the low range 
of temperature^^Ti^^. In presence of magnetic field, the 
oscillation of the diffusive thermopower has been studied 
theoretically as well as experimentally^""—. It is seen in 
the low magnetic field regime that both Sxx and Sxy are 
periodic in inverse of the magnetic field. This is due to 
the oscillating density of states of the 2DEG in presence 
of magnetic field. 

There is no theoretical or experimental study on mag- 
netothermoelectric properties of the 2DEG systems with 
the Rashba SOI. We report here for the first time the 
effect of the Rashba SOI on thermal transport proper- 
ties of a 2DEG in presence of perpendicular magnetic 
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field. The total thermalconductivity and the total ther- 
mopower produce beating patterns because the thermo- 
electric coefficients for spin-up and spin-down electrons 
oscillate with two closely spaced different frequencies. By 
analyzing the beating pattern, we find a simple equa- 
tion which determines the Rashba SOI strength if the 
magnetic fields corresponding to any two successive beat 
nodes and the number of oscillations in between are 
known from the experiment. 

This paper is organized as follows. In section II, we 
briefly mention the energy spectrum and the DOS of the 
2DEG with the Rashba SOI for zero and non-zero mag- 
netic field cases. In section III, we have studied the ther- 
moelectric coefficients for zero magnetic field case. We 
also provide the formalism to be used for studying ther- 
moelectric coefficients in presence of magnetic field. In 
section IV, we present our numerical and analytical re- 
sults. We provide a summary and conclusion of our work 
in section V. 



II. ENERGY SPECTRUM AND DENSITY OF 
STATES OF A 2DEG WITH THE RASHBA SOI 

A. Zero magnetic field case 

The Hamiltonian of an electron with the Rashba SOI 
is given by^ 

^2 

(1) 



where p is the two-dimensional momentum operator, m* 
is the effective mass of the electron, 1 is the unit matrix, 
c = (ca;, cTy, Cz) are the Pauli spin matrices and a is the 
strength of the Rashba SOI. At non-zero momentmn, the 
spin degeneracy is lifted due to the presence of the SOI. 
The energy spectrum of the "spin- up" and "spin-down" 
electron is given by 



2m* 



±a\k 



(2) 



Here, the -I- and - signs correspond to the spin-up and 
spin-down electrons. The density of states (DOS)^ for 
spin-up and spin-down electrons are 
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Here, Do = m* /{-kK^), E^ — 2m*o? jf? is the Rashba 
energy determined by the Rashba SOI strength a and 
0(i?) is the unit step function. 



B. Non-zero magnetic field case 

The Hamiltonian of an electron (— e) with the Rashba 
SOI in presence of a perpendicular magnetic field B = Bz 
is given by 
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^g^^Ba^, (5) 



where /i^ — eh/{2me) is the Bohr magneton with rrie 
is the free electron mass and g is the effective Lande g- 
factor. The exact energy spectrum and the corresponding 
eigenfunctions of the above Hamiltonian are derived in 
Ref.^^. The resulting eigenstates are labeled by a new 
quantum number s. For s = 0, there is only one energy 
level which is same as the lowest Landau level without the 
Rashba SOI. The corresponding energy is given by E^ = 
Eq — {hu! — giigB)/2. Here, lo = eB/m* is cyclotron 
frequency. For s = 1,2,3...., there are two branches of the 
energy levels, denoted by + corresponding to the "spin- 
up" electrons and — corresponding to the "spin-down" 
electrons with energies 



Ef = shhj±\ El + sEafujj. 



(6) 



Using the Green's function method, the DOS for spin- 
up and spin-down electrons in presence of magnetic field 
are calculated in Ref.— . These are given by 
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where Fq is the impurity induced Landau level broaden- 
ing. 

III. THERMOELECTRIC COEFFICIENTS 

In this section, we shall develop the formalism for the 
thermoelectric coefficients of a 2DEG with the Rashba 
SOI system for both the cases: zero and non-zero mag- 
netic fields. 



A. Zero magnetic field case 

In this sub-section, we consider a 2DEG with the 
Rashba SOI and calculate the thermal power and thermal 
conductivity. Within the linear response regime, the elec- 
trical current density J and the thermal current density 
J„ for spin-up and spin-down electrons can be written as 



and 



(8) 
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where E is the electric field and L!j! with i,j = 1,2 are the 
phenomenological transport coefficients for spin-up and 
spin-down electrons in absence of magnetic field. These 
are the main equations that determine the response of 
a system to the external forces such as electric field and 
temperature gradient. In presence of the Rashba SOI, 
the spin-up and spin-down electrons will contribute to 
the total electrical and thermal current. Therefore, the 
total electrical current and the thermal current densities 
are 



We calculate the total electrical conductivity <j{Ep) at 
the Fermi level, which is given as 



where (Tq 
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* 9 9 

TO e To a 



1 



(17) 



is the Drude conductivity without SOI. 
The similar expression of the Drude conductivity is ob- 
tained by using a different method in Ref.^^. The total 
thermal conductivity is then 



and 



(10) 
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Here, = Lf + L^_l and L^^ can be written in terms 
of the integral /M; = /W,^^! = tL^^ = -I^^^e, 
= /(2)/(e2r). Also, /W = -K/M^- with 
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where r = 0, 1, 2 and f{E) = 1/[1 + exp{E - ri)l3] is the 
Fermi-Dirac distribution function with rj is the chemical 
potential and (3 = l/(fc^T). Here, cr~^{E) and cr~{E) are 
the energy-dependent conductivity for spin-up and spin- 
down electrons, respectively. In an open circuit condition 
(J = 0), the thermopower is given hy S = L^"^ / L^^ . 
Then at low temperature, diffusion thermopower S and 
the diffusion thermal conductivity k can be expressed in 
terms of the electrical conductivity through the Mott's 
relation and the Wiedemann-Franz law as 



dE 



\na{E) 
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and 
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We note that the thermal conductivity and the ther- 
mopower is enhanced due to the presence of the Rashba 
SOI. 



B. Non-zero magnetic field case 

In this subsection, we shall study the thermoelectric 
coefficients of a 2DEG with the Rashba SOI in presence 
of the perpendicular magnetic field. Thermoelectric coef- 
ficients in presence of magnetic field (without SOI) were 
obtained by modifying the Kubo formula in ReL— i— . 
Here we shall generalize these results to the SOI sys- 
tems. These phenomenological transport coefficients can 
be re-written as 



/•(o),± 
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<^ = ^[^l^>^ - eTiC'^'^^^S^),.], (21) 
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Here, Lq — (vr^fc^ )/(3e2) is the Lorentz number and 
a{Ep) = a^{Ep) + a~{Ep) is the total electrical con- 
ductivity at the Fermi energy. 

By using the Boltzmann transport equation, we eval- 
uate the zero-temperature energy-dependent electrical 
conductivity for spin-up and spin-down electrons, which 
are given by 



a^{E) = —r{E)g^{E) 



(15) 



Assuming the energy dependent scattering time to be 
T — To{E/Ep)P , where p is a constant depending on the 
scattering mechanism. We also assumed that r is the 
same for spin-up and spin-down electrons. Substituting 
Eqs. ©, O and ^ into Eq. i^^, then the diffusion 
thermopower is obtained as 
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where 



= / dE 



df{E) 
dE 



{E~vYai^{E). (22) 



Here, n^v = x,y. Also, (j^jj{E), S^^, and are the 
zero-temperature energy-dependent conductivity, ther- 
mopower and thermal conductivity tensors, respectively, 
for spin-up and spin-down electrons. The total ther- 
mopower and thermal conductivity can be obtained from 
Sfii, = 5*+^ + S'^jj^ and = + k^^. 

In electron systems, conduction of carriers takes place 
by the diffusive and collisional mechanisms. The colli- 
sional contribution leads to the SdH oscillation with in- 
verse magnetic field due to the quantized nature of the 
energy spectrum. We will consider the collisional mech- 
anism only because electrons do not possess any drift 
velocity in our case. In the linear response regime, the 
conductivity tensor can be written as the sum of di- 
agonal and non-diagonal as a^i, = a"^^ 
(t"^^ is the Hall contribution. Here, a-x- 



(t"^ , where 
= cr™' and 
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(T™ . Similarly, for the thermal trans- 
port coefhcients the following relations are valid: Cxx = 

£(r)col ^(r)col i fix) /•(^}clif ^(r)col ^(r)col 

XX ~ ^yy anct l~'yy — *^yy ~r ^yy — ^yy 

The exact form of the finite temperature collisional con- 
ductivity has been calculated in RefJ^ for the screened 

impurity potential C/(q) = 27re^/ + 9y + ^s) in mo- 
mentum space. Here, ks is the inverse screening length 
and e is the dielectric constant of the material. In the 
limit of small |q| <C ks, C/(q) — 27re^/(efcs) = Uq. In this 
limit, one can use Tq ~ nPh^ /NjUq with tq is the colli- 
sional time, / = y^h/eB is the magnetic length scale, Uq 
is the strength of the screened impurity potential and Nj 
is the two-dimensional impurity density. The exact form 
of the finite temperature conductivity can be reduced 
to the zero-temperature energy-dependent electrical con- 
ductivity as 



(23) 



h 2ttToP ' 

where if = [(2s T - 2sDl + (2s ± l)]/^^ with 

Ds = y^sEahuj/lEo + + sEJ}bj\ and As = \^D\ 



Using Eq. ([22|) . the finite temperature diagonal and off- 
diagonal coefficients (£2 and c}yx) can be written as 



£(r),± ^ 



and 



E=Ef 

(24) 



(25) 



IV. 



NUMERICAL RESULTS AND 
DISCUSSIONS 



In our numerical calculations, the following parame- 
ters are used: carrier concentration rig = 3 x 10^^ /m^, 
effective mass m* = O.OStoe with me is the free electron 
mass, 5 = 4 and the Rashba SOI strength a — 5 x 10~^^ 
eV-m and To = 0.01 meV. For better visualization of the 
oscillations, we have used T = 1 K for the thermopower, 
T — 0.5 K for the thermal conductivity. In Fig. [1], 
the components of the thermopower tensor in units of 
—kgje are shown as a function of the inverse magnetic 
field. The diagonal thermopower components Sxx and 
Syy are identical and therefore only Sxx is shown. In 
Fig. [2], the thermal conductivity is shown as a function 
of the inverse magnetic field. The magnetic field depen- 
dence of the thermal conductivity is same as that of the 
electrical conductivity. Figures [1] and [2] show the ap- 
pearance of the beating pattern in the thermopower and 
thermal conductivity. 




FIG. 1: Plots of the thermopower versus inverse magnetic 
field. In the lower panel, dashed and soUd lines correspond to 
the analytical and exact results, respectively. 




1/B{1/T) 

FIG. 2: Plots of the longitudinal component of thermalcon- 
ductiviy versus inverse magnetic field B. The dashed and 
solid lines correspond to the exact numerical and analytical 
results. 



To analyze the beating pattern in the thermoelectric 
coefficients, we shall derive analytical expressions of the 
thermoelectric coefficients. The components of the ther- 
mopower for spin-up and spin-down electrons are given 
by 



S: 
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XX 
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The dominating term in the above two equations is the 
last term. The analytical form of k^^ and S^y can be 
obtained directly by deriving analytical form of the phe- 
nomenological transport coefficients. The analytical form 
of the DOS given in Eq. ^ allows us to obtain asymp- 
totic expressions of Sxy and Kxx- This is done by replac- 
ing the summation over discrete quantum numbers s by 
the integration i.e; — )■ 27rP J D^{E)dE, then we get 



and 



/:(2),± ^ (I 



8(wro)2 



1 _ ^ 
3 2 



G"(a;)cos (Ztt^) 



where the impurity induced damping factor is 



(29) 



(30) 



and the temperature dependent damping factor is the 
derivative of the function G{x) with G{x) = xj sinh(a;). 
Here, x = T/Tc with Tc = huj/2iT^k^. Note that Gix) is 
the temperature dependent damping factor for the elec- 
trical conductivity tensor. Also, the oscillation frequen- 



cies are 



he 



Ef + ^T ^Je^+E^Ef 



(31) 



The off-diagonal thermopower Sxy for spin-up and 
spin-down electron is obtained as 



9* = 



^dG'{x) sin(27r^) 



e 4ajTo 

The total thermopower is given as 

Sxy = -^-^noG'ix) sin (2^4) cos (2n^ 
e 2uJTo \ B J \ B 



(32) 



(33) 



Here, U = (/+ + D/^ and U = (/+ " /")/2. ^ In 
the lower panel of Fig. [1], we compare the analytical 
expression of Sxy with that of the numerical result. The 
analytical result matches very well with the numerical 
results. 

For thermal conductivity, the dominant term 
is C^xx- The approximate analytical form of can be 
obtained from Eqs. (PT|) and ((^ as 



Ln 



Tao 



1 - ^nDG"{x) cos f27r4- 
2 \ B 
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The total thermal conductivity can be written as 



(34) 
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(35) 



Equations (|32)) and ([34| show that the thermopower 
and the thermalconductivity of spin-up and spin-down 
electron oscillates with different frequency and /~, 
respectively. Therefore, the beating pattern appears in 
the total Sxy and n. It is quite difficult to obtain the 
analytical expression of Syy, but the origin of the oscil- 
latory part is due to the oscillatory density of states at 
the Fermi energy. 

We get the condition for beating nodes from 
the periodic term with frequency difference /d: 
cos{2tt fi/ B) B=Bj — which gives 



A^, + {l~g*Y{hu:,Y = nw,{j + ]-). (36) 



Here, = 2k is the zero-field spin splitting energy 
with is the Fermi wave vector, j — 1, 2, 3.. is the j-th 
beat node and g* = gm* /{2me). Also, ujj — eBj/m* and 
Bj is the magnetic field corresponding to the j-th beat 
node. Using the above equation, one can determine the 
zero-field spin splitting energy or the Rashba strength if 
we know the number (j) of any node and the correspond- 
ing magnetic field Bj. In practice, the numbering of the 
beat nodes is quite difficult. The above equation can be 
re-written for two successive beating nodes as 



(37) 

Therefore, the Rashba SOI strength can be determined 
from Eq. (|37p by knowing the magnetic fields correspond 
to any two successive beat nodes. 

In the above analytical expressions [Eqs. ([55)1 and ([55]) ] 
the periodic term with frequency /a gives the number 
of oscillations between the two successive beat nodes as 
given by 



E^ 



1 



B 



Therefore, we can also determine the Rashba strength 
from Eq. (j38l) by knowing the magnetic fields correspond 
to any two successive beat nodes and the number of oscil- 
lations in between. We note that Eqs. ((36|) and ([38| are 
the same as obtained in the beating pattern formation in 
the SdH oscillations^^. 



V. CONCLUSION 

We present theoretical study of the effect of the Rashba 
SOI on the thermoelectric coefficients. In absence of mag- 
netic field, the thermopower and the thermal conductiv- 
ity are enhanced due to the presence of the SOI. The nu- 
merical results of all the thermoelectric coefficients are 
given. In addition to the numerical results, we provide 
the analytical expressions of the off-diagonal component 
of the thermopower (Sxy) and the diagonal components 
of the thermal conductivity (kxx)- The appearance of 
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the beating pattern in the thermoelectric coefficients can 
be explained from the fact that the two branches oscil- 
late with slightly different frequency and produce beating 
pattern in the thermoelectric coefficients. The analytical 
results match very well with the numerical results. The 
Rashba SOI strength can be determined if the magnetic 
field corresponding to any two successive beat nodes are 
known from the experiment. 
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